r^-Exponents in the one dimensional antiferromagnetic Heisenberg model with next to 

nearest neighbour coupling 

C. Gerhardt, A. Fledderjohann, E. Aysal, K.-H. Mutter* 
Physics Department, University of Wuppertal, 4^097 Wuppertal, Germany 

J.F. Audet, H. Kroger 

Department of Physics, Universite Laval, Quebec, Canada 

(February 1, 2008) 

We investigate the critical exponents 7)3(0, M), 7)i(a, M) associated with the singularities in the 

longitudinal and transverse structure factors of the one dimensional antiferromagnetic Heisenberg 

model with nearest (Ji) and next to nearest (J2) neighbour coupling of relative strength a = -^ 

^\ ' and an external field B with magnetization M{B). 
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I. INTRODUCTION 

r-j \ In this paper we continue the investigation of the one dimensional spin i Heisenberg model 

in 

H = 2^{S{x)S{x + 1) + aS{x)S{x + 2)) + 2B ^ Saix) (1.1) 

'""' 1 x=l x=l 

> 

CO ' with next to nearest neighbour coupling parameter a and external field B. 

00 . Let us briefly summarize those results [ |l| - [ H relevant for our later investigation. 

(N 

(^ In the absence of an external field {B = 0), the ground state of the model is a singlet {S = 0) state with momentum 
^\ pq, where po = for iV = 4,8, 12, ... and po = tt ioi N = 6, 10, 14, .... Tliia statement holds at least in the interval 
On ~k "^ '^ ^ 1- AtQ; = ^ the model reduces to the Majumdar Ghosh niodelE2l with degenerate dimer ground states, 
^i The quantum numbers of the first excited state change with a: There is a triplet S' = 1 state for < a < ac = 0.241... 
2 '. and a singlet (5 = 0) state for ac < a < ^. The momentum of the first excited state is pi = po + tt. Moreover there 
is no gap in the 'spinfluid' phase and a gap in the 'dimer' phase a > Uc- The structure of these two phases can be 
exploited by means of the static and dynamical correlation functions of appropriate operators. In the spinfluid phase 
< a < ac the AS* — 1 operator 



X3 
O . 1 

O. ^3(p)^ ^e^p-53(:.) (1.2) 



^^ ■ generates the transition from the singlet to the triplet excited state. To examine the dimer phase a > ac we need a 
Vh ■ AS* = operator. These transitions are generated by the dimer operator 



Dip) = ^J2 e'''^0ix)Six + 1)- < Six)S{x + 1) >) (1.3) 



The corresponding static structure factors 

S3{a,p,N) =< S+{P)S3{P) > and D{a,p,N) =< D+{p)D{p) > (1.4) 

behave as follows for A^ — > 00: 

S3{a,p = TT, A^) diverges logarithmically for a < ac but stays finite for a > ac- 

D{a,p — TT, N) diverges with a power depending on a for a > ac and stays finite for a < ac- The power behaviour 

degenerates to a logarithmic behaviour for a — ac- 

In the presence of an external field B the ground state of the model has total spin S = M ■ N, where M = M{a, B) 
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is the rnagnetization. The behaviour of the magnetization curve M{a, B) near saturation B ^ Bg, M ^ ^ changes 
with aH. J-, 

At a = it is known to develop a square root singularityll3: 

Mia^O,B)^^--iB,-B)K iov B ^ B„ (1.5) 

whereas the numerical data for a — j{N < 28) support a quartic root singularity for a = ^ M: 

Mia=],B)^l-^iB,-B)i £4 = 1.70(5), for B ^ B,. (1.6) 

The gap in the dimer phase a > ac appears in the low field behaviour of the magnetization curve: 

M{a, B)^0, B < Bc{a), a > ac- (1.7) 

Note that the model with B > Bc{a) is gapless, provided that there are no 'plateaus' {M(a,B) = const, Bic < B < 
B2c) in the magnetization curve. l— ■ 

The momentum of the ground state ps follows Marshall's sign rulal3: 

ps = for 2S' + iV = 4n, ps = tt for 25" + iV == 4n + 2 (1.8) 

for < a < i. 

The singularities in the static structure factors change, if we switch on an external field: The transverse structure 

factor at p = TT, a = 

Si{a,p = TT,M,N)KiBi{a,M)N^-'^''^°''^'^^ +Ai{a,M) (1.9) 

diverges with a field dependent critical exponent 771 (a, M). ?7i(a = 0,M) has been calculated in Ref. 04 by means 
of the Bethe Ansatz. A second, weaker singularity, which moves with the external field appears at the softmode 
momentum p = pi{M) — 2'kM. To our knowledge, the positions of both singularities do not depend on a. The 
critical exponent r]i{M,a) however does. It has been found to beE£l 

77i(a = 0,M= -) =0.65 and 77i(a = -, M = -) = 1.16 (1.10) 

The longitudinal structure factor at p = tt (Af > 0) stays finite for a < | but develops a singularity at the field 
dependent softmode momentum p = P3{M) = 7r(l — 2Af ): 

S3ia,p = P3(M), M, N) « Bsia, M)N'-'''<-"'^'^ + A^ia, M) (1.11) 

Again, the position of the singularity does not depend on a, whereas the critical exponent rj^ia, M) changes drastically 
with a: 

rjsia = 0, M = t) = 1-50 and 7/3(0 = -, Af = -) = 0.84 (1.12) 



The soft mode singularity at p = p^jAI ) is a lso found in the dimer structure factor D{a,pJ^) defined in (1.3) and 
(1.4). A finite size analysis of the type ( |l . 1 1 ) for D{a,p, N) yields for the critical exponentalj: 



rjoia = 0, M = ^) = IA9 77^(0 = i, M = i) = 0.82 (1.13) 



These values almost coincide with those of the longitudinal structure factor given in ( 1.12| ) 



It is the purpose of this paper to determine the complete a dependence of the critical exponents 771(0;, M), 773(0;, Af). 
The paper is organized as follows: 

In section 2 we exploit the range of validity of Marshall's sign rule ( |1.8| ). In section 3 we study the impact of the 
next to nearest neighbour coupling on the lowlying excitations and on the static structure factors. The latter are 
computed numerically on systems up to iV = 32. A finite size analysis of ( |l.9| ) and ( |l.ll[ ) yields the critical exponents 
771(0, Af), 773(0, Af). Section 4 is devoted to the study of an unexpected phenomenon, which we found for negative 
o- values: For o < a-{M) the finite size behaviour of the longitudinal structure factor ( jl.ll ) changes systematically 
from a monotonic increase to a decrease. 



II. THE IMPACT OF FRUSTRATION ON THE GROUND STATES. LEVELCROSSINGS. 

It was pointed out in the introduction that the momenta psia) of the ground state \S,Sz — S,ps{a) > in the 
sectors with total spin S foUow MarshaU's sign rule for a < j. We found deviations from this rule for 

a > a,{M - |) (2.1) 

We computed the ground state energies E{a,ps{a),M = j^,N) on small systems with N — 10, ..,20 sites. The 
dependence on the frustration parameter a is shown in Fig. 1 for N = 12. Here we have marked the different ground 
state momenta by different symbols. 
At Af = -^ = deviations from Marshall's sign rule occur first at: 

ao{M = 0) = i (2.2) 

Here we meet the Majumdar Goshll3 model, which is known to have two degenerate ground states, namely dimer 
states with momenta p — and p — n, respectively. 

A twofold degeneracy - with respect to the ground state momenta pg {a),pg (a) emerges along the whole curve 
a = ao{M), which is plotted for N=10,12, 14,16,18,20 in Fig.2 



.(1)^ 



The first momentum pg (ao) follows Marshall's sign rule (L8). We have looked for an empirical rule for the second 

momentum Pg {cto) but we did not find such a rule which holds for all momenta and system sizes N. 
In the saturating field limit M -^ i the curve ao(M) meets the point 

MM -^l) = \ (2.3) 

Indeed, the eigenvalue problem can be solved analytically for 5* = — — 1 with the ansatz (1 magnon states) 



N 

where \x > denotes a spin state with spin — ^ at site x and spin ^ at all other sites. The energy of the 1 magnon 
state is found to be: 

11 N 

E{a,p, M=---,N)^2 cosp + a2 cos(2p) + (_ - 2)(H- a) (2.5) 



The ground state energy and its momentum p — ps{a) follows by minimizing (2^) with respect to p. For a < j 
the ground state momentum is found to be pg (a) = tt, 5* = -y — 1 in accord with Marshall's sign rule. For a > j 
however, the minimum is found for p = Pg (a) where 

cosp^^Ua)^— S = 1 N^oo (2.6) 

On finite lattices, the difference between the two momenta turns out to be: 

AP5(«0 = \) = bi'H^o = J) -Pf («o = J)l = ^ ^ = Y " ^ (2-^) 

As a consequence of the levelcrossing at a = ao(A/ — j^), M fixed, the derivatives of the ground state energies 

^E{a,ps{a),M = ^,N) (2.8) 

change discontinuously, as can be seen in an amplification of Fig. 1. 



III. SOFTMODES IN THE EXCITATION SPECTRUM AND THE ASSOCIATED J7-EXPONENTS. 

We have studied the finite-size dependence of the energy gaps 

S + AS S 

WAs(a,P, M, N) = E{a,p = ps + P, M = ^ , N) - E{a,p = ps, M = -,N) (3.1) 

for AS = and AS" = 1 in the domain a < ao{M) where the ground state momentum follows Marshall's sign rule. 
In this regime the gap ti-'AS=i(a,P = tt, M, iV) vanishes in the thermodynamical limit in such a way that the scaled 
quantity 

lim 7VtJAS=i(a,P = 7r,M,7V) = f7i(a,M) a < ao(M) (3.2) 

converges to a finite non- vanishing limit. The same holds for the gap was=o(q;,p = P3(M), M, N), pj,{M) = 7r(l — 2Af ), 
if the next to nearest neighbour coupling a is positive: 

lim NujAS=oia,P^P3iM),M,N)^n3{a,M) (3.3) 

For negative a-values (a < a^{M) < 0) however, we observe a tendency in the numerical data, which at least hints 
to the emergence of a gap at the momentum p — p3(M): 

lim WA5=o(a,P = P3(M),Af,7V) = A3(a,M) a < a_(Af) < (3.4) 

as can be seen from Fig. 3 . It is hard to decide from the finite system results (A^ = 16, 20, 24) the exact position 



a = a-{M) where the gap (3.4) opens 



In the gapless regimes, where ( |3.2[ ) and ( |3.3| ) is valid, we expect that the critical behaviour of the system is properly 
described by conformal field theory. This means in particular that the ratios [ |l^ Q 

20,(a,M) = ^M^, a = 3,1, (3.5) 

7rw(a, M) 

can be identified with the critical exponents r]a{a, M): 

2ea{a,M)^Tja{a,M), a = 1,3. (3.6) 

Here 

z;(a,M) = -^lim iV(iJ(a,p = ps + ^, Af = ^, TV) - i?(a,ps, A/ = ^, ^)) (3.7) 

is the spinwave velocity. For fixed valu es of M {M = i i , i ) we have determined the a-dependence of 26a (a, Af ) , a ■ 



1,3, from the energy differences (3J_) (3/7) as they enter in the ratios (3^). The result can be seen from the solid 



curves in Figs 4a,b,c. The solid dots repr esent the determination of the critical exponents 771 (a. A/), 773(0, M) as they 



follow from a fit of the form (L£), ( I.ll[ ) to the finite system results {N < 32). Comparing the two determinations 



we come to the following conclusions: 

1. The identity 29i{a, M) = 77i(a, M) for the critical exponent in the transverse structure factor is well established 
for -0.5 <a< 0.25. 

The same holds for the identity 26*3(0;, Af) = 773(0, Af) for the critical exponent in the longitudinal structure 
factor in the interval a-{M) < a < 0.25. If we approach the curve o = oo(Af ) the convergence of the Lanczos 
algorithm slows down more and more, due to the emergence of the level-crossing discussed in section 2. 

2. The two curves 20i(o, Af ), 26*3(0, Af ) cross each other at o = Oc(Af ) 

20i(oe(Af), Af) = 203(ac(Af),Af) = 20(Af) (3.8) 

where 

a,(A/ = i) = 0.18 Oe(Af = J) = 0.20 Oc(Af = ^) = 0.24 (3.9) 



and 

26i(Af = -) = 1.01 20(Af = -) = 1.02 26'(Af = -) = 1.02 (3.10) 

The a-values are quite close to the transition point a^M = 0) = 0.241 from the spinfluid to the dinier phase. 
The same holds for the critical exponents ri{M), which deviate only slightly from ri{AI = 0) = 1. 

3. The relation 

40i{a, M)e3{a, M) = 1 (3.11) 

appears to be satisfied within a few percent for < a < ^ 

4. For negative values of a, we observe in the data for 773 (a, M) a discontinuou s stru cture (open symbols). Looking 



at the numerical data, which enter in the determination of 773(0, M) via eq. (1.11), we found a systematic change 
in the finite size dependence. For a > a^{M) 

a_(Af = i) =-0.31 a_(Af = i) = -0.19 a_(Af = i) = -0.15 (3.12) 

the longitudinal structure factor monotonically increases with N, whereas it decreases for a < a_(Af). In the 



latter regime we expect the emergence of the gap (3.4) 



IV. THE DISAPPEARANCE OF A FIELD DEPENDENT SOFTMODE 



The change in the finite size dependence of the gap (3.3, 3.4) and of S3{a,p3{M), M, N) provides us with a first 



hint, that the field dependent softmode at p = PsiM) = 7r(l — 2Af) might disappear for a < a_{M). In this section, 
we are looking for further evidence for this hypothesis. In Figs 5a, b we compare the momentum distribution of 
S3{a,p,M = i,iV) for a = a-{M = |) = -0.19 and a = -0.4, respectively. At a = a-{M = j) = -0.19 (Fig 
5a) the momentum distribution is well approximated by two straight lines with different slopes for p < 733 (Af) and 
p > ps^M), respectively. This discontinuity is more and more washed out, if the next to nearest neighbour coupling 
decreases further. E.g. at a = —0.4 (Fig 5b) the p-distribution of the longitudinal structure factor appears to be 
smooth in the thermodynamic limit. The approach to this limit is indicated by an arrow. 83(0, p, M, N), a < a-{M) 
is monotonically decreasing with N for p < p^lM) but increasing for p > p3{M). 
A more drastic effect can be seen in the dynamical structure factor: 

S3{.a,Lo,p,M,N)=Y,5{LO-{En-E,))\<n\S3{p)\s>\^ (4.1) 

n 

which we computed by means of the recursion methodcj cA for M — j and A^ = 28. The excitation spectrum is 
plotted in Figs 6 a,b for a = —0.19 and a = —0.4, respectively. The numbers denote the corresponding relative 
spectral weight in percentage terms. The curves guide the eye to the excitations with dominant spectral weight. For 
a = —0.19 (Fig 5 a) the spectral weight is distributed over a band of excitation energies which broadens in the vicinity 
of the momentum p = 733 (Af). For a = —0.4 (Fig 5b), however, the spectral weight is more concentrated at higher 
excitation energies. In particular, the lowest excitation at p = P3{M) has a relative spectral weight less than 10% for 
AT = 28. 

V. DISCUSSION AND CONCLUSION 

In this paper, we studied the impact of a next to nearest neighbour coupling a and an external field B on the zero 
temperature properties of the one dimensional spin i antiferromagnetic Heisenberg model. We found the following 
features: 



1. The momentum of the ground state follows Marshall's sign rule ( |1.8| ) for a < ao{M) where ao(0) = ^ and 
(20(2) = 4- The ground state is twofold degenerate with respect to its momentum for a = ao{M) 



A study of the finite size dependence (1.9) and (1.11) yields the a-dependence of the critical exponents 



rii{a, M),r]3{a, M) associated with the softmode singularities at p = tt and p = p^^M) i n th e tran sverse and 



longitudinal structure factor, respectively. Good agreement is found with the prediction (3.5), (p.q) of confor 



mal field theory for r]i{a,M) — ^ < a < j and for 773(0, M) < a < j (Figs 4a-c). For these g-v alues 



the spectral weight - entering into the definition of the corresponding dynamical structure factors (4.1) - is 
concentrated around the lower bound of the excitation spectrum. This seems to be a crucial condition in or- 
der that the critical behaviour is described correctly by conformal field theory. In the thermodynamical limit 
the dynamical structure factors Si{a,uj,p = tt, M) and 83(0, to, p = p3{M), M) develop infrared singularities 
ijj-('^-Va(a)) a = 1^ 3^ which can clearly be seen in a finite size scaling analysis. Such an analysis was performed 



in Ref. 14 for a = 0. 



3. Deviations from the relation (3.7) - predicted by conformal field theory - appear in the longitudinal case a = 3 



(Figs 4a-c) for negative values of the next to nearest neighbour coupling and increasing Af-values. This is 



accompanied by the fact that the spectral weight in (4.1) is distributed over a band of excitation energies, which 
broadens with decreasing values of a. 

4. There are several indications that the field dependent soft mode at p = P3{M) = 7r(l — 2M) disappears for 
negati ve nex t to nearest neighbour couplings a < a_(M) < 0: A gap (BA) opens and the longitudinal structure 



factor ( 1.11 ) changes its finite size dependence from a monotonic increase to a decrease. Moreover the cusp-like 
singularity in the momentum dependence at p = p3(M ) is washed out and the spectral weight is shifted from 
low to higher excitation energies. 

Therefore, we find a further confirmation of the hypothesis - formulated in Ref. G2 - namely that field dependent 
soft modes only exist if the system is sufficiently frustrated. As was pointed out in Ref. ^ this condition is not 
satisfied in the twodimensional spin ^ antiferromagnetic Heisenberg model with nearest neighbour coupling. 
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Figure Captions 

FIG. 1. The ground state energies E{a,ps,M — jf,N) in the sector with total spin S on a ring with A'' = 12 sites. The 
ground state momenta ps(q) change with the next to nearest neighbour coupling a as indicated by the different symbols. 

FIG. 2. The curve a = ao{M = jf) where the ground state |S, p > in the sector with total spin S is degenerate with respect 
to the momentum p = pg {a),pg (a). 

FIG. 3. Finite size dependence of the gap a;A5=i(Q,P3(M), M,iV) for a = 0.1, 0.0, -0.1, -0.2, -0.3, -0.4, -0.5. 



FIG. 4. Comparison of t he ratio 2 6i(a, M) i = 1, 3 ( |3.5| ) (solid curves) and the critical exponents 77i(a, M) i = 1, 3 in 
the static structure factors (p^( |l.ll| ). a.) M = \ h) M = \ c) M ^ \ 

FIG. 5. The momentum dependence of the longitudinal structure factor Sz[a,p,M — \,N) N = 28,24,20,... . a) 
a = a_(M = i) = -0.19 b) a = -0.40 

FIG. 6. Excitation energies and relative spectral weights in the dynamical structure factor Sz{a,uj,p,M = j,N = 28) in 
percentage terms. The lines connect the excitations with the dominant spectral weight, a.) a — a-{M = j) = —0.19 b) 
a = -0.4 
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